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Abstract 

A  three-dimensional  numerical  model  is  developed  to  study  the  electrochemical  species  characteristics  in  a  free-breathing  cathode  of  a  direct 
methanol  fuel  cell  (DMFC).  A  perforated  current  collector  is  attached  to  the  porous  cathode  that  breathes  the  fresh  air  through  an  array  of  orifices. 
The  radius  of  the  orifice  is  varied  to  examine  its  effect  on  the  electrochemical  performance.  Gas  flow  in  the  porous  cathode  is  governed  by  the 
Darcy  equation  with  constant  porosity  and  permeability.  The  multi-species  diffusive  transports  in  the  porous  cathode  are  described  using  the 
Stefan-Maxwell  equation.  Electrochemical  reaction  on  the  surfaces  of  the  porous  matrices  is  depicted  via  the  Butler- Yolmer  equation.  The  charge 
transports  in  the  porous  matrices  are  dealt  with  by  Ohm’s  law.  The  coupled  equations  are  solved  by  a  finite-element-based  CFD  technique.  Detailed 
distributions  of  electrochemical  species  characteristics  such  as  flow  velocities,  species  mass  fractions,  species  fluxes,  and  current  densities  are 
presented.  The  optimal  breathing-hole  radius  is  derived  from  the  current  drawn  out  of  the  porous  cathode  under  a  fixed  overpotential. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Passive  direct  methanol  fuel  cells  (DMFCs)  have  received  a 
great  deal  of  attention  since  they  are  considered  to  be  a  poten¬ 
tial  candidate  in  replacing  batteries  in  portable  electronic  device 
applications  [1-5].  In  structure,  they  are  much  simpler  than 
active  ones  (which  have  fuel  pumps  and  oxidant  supplies)  and, 
therefore,  the  parasitic  power  penalty  from  the  auxiliary  devices 
is  eliminated.  However,  the  power  density  from  passive  systems 
is  not  as  good  as  those  running  in  an  active  mode  due  to  poor 
transport  [6-11].  It  is  well  known  that  an  insufficient  oxygen 
supply  to  the  active  sites  will  result  in  uneven  current  distribu¬ 
tions  in  the  electrode  which  reduces  the  fuel  cell  performance. 
Thus,  one  of  the  challenges  in  the  design  of  passive  DMFCs  is 
how  to  feed  the  electrode  with  reactants  and  remove  the  products 
from  the  electrode  appropriately. 

In  this  paper,  the  effect  of  the  breathing-hole  size  on  the  cath¬ 
ode  performance  of  a  passive  DMFC  is  studied.  A  staggered 
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arrangement  of  circular  holes  is  distributed  on  the  current  col¬ 
lector  that  is  attached  to  the  porous  cathode.  Fresh  air  passes 
into  the  porous  cathode  through  these  holes  while  the  electrical 
current  is  conducted  out  through  the  solid  counterpart.  Conserva¬ 
tion  equations  of  coupled  mass,  momentum,  species  and  charge 
are  developed.  The  momentum  conservation  in  the  porous  cath¬ 
ode  is  described  by  Darcy’s  law.  The  multi-species  diffusion  in 
the  porous  cathode  is  illustrated  by  the  Stefan-Maxwell  equa¬ 
tion.  The  Butler- Volmer  equation  describes  the  electrochemical 
reaction  on  the  surfaces  of  the  porous  cathode.  A  finite-element 
based  computational  fluid  dynamics  methodology  [12,13]  is 
employed  to  solve  the  multi-physics  transports  for  the  free- 
breathing  cathode  of  a  DMFC.  Detailed  species/electrochemical 
characteristics  such  as  the  flow  velocity,  species  mass  flux, 
species  mass  fraction,  current  density,  and  overpotential  are  pre¬ 
sented.  The  electrochemical  performance  of  various  perforation 
sizes  is  derived  from  the  amount  of  current  that  is  drawn  out 
of  the  porous  cathode  under  a  fixed  total  overpotential.  The 
results  obtained  in  the  present  study  can  help  in  understand¬ 
ing  the  local  gas  transport  and  electrochemical  characteristics 
of  a  free-breathing  DMFC.  Also,  they  provide  a  solid  basis  for 
optimizing  the  geometry  of  the  passive  DMFC  stack. 
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Nomenclature 


Ci 


Dh2o 

Dq2 

F 

h  t 
M 

n 

P 

Pi,m 


Pi,s 


r 

R 

T 

u 

x9  y,  z 


mole  concentration  of  the  species  i  (molm  3), 


p 

RT 


diffusivity  of  water  vapor  (m2  s-1) 
diffusivity  of  oxygen  (m2  s-1) 

Faraday’s  constant 
current  density  (A  m-2) 
molecular  weight  (kg  mol-1) 
mass  flux  (kgm2  s-1) 
pressure  (Pa) 

possibility  of  the  electrolyte  in  the  connection  of 
the  catalyst  layer 

possibility  of  the  catalyst  in  the  connection  of  the 

catalyst  layer 

radius  of  the  perforation 

universal  gas  constant  (W  mol-1  K-1) 

temperature  (K) 

velocity  vectors  (ms-1) 

coordinate  system,  Fig.  1  (m) 


Greek  symbols 
a  symmetric  factor 

e  cathode  porosity 

rjt  total  overpotential  across  the  computational 
module  (V) 

k  permeability  (m2) 

ft  dynamic  viscosity  (m  s-2) 

p  density  (kgm3) 

o' i  ionic  conductivity  of  the  membrane  phase 

(S  m-1) 

crs  electric  conductivity  of  the  catalyst  phase  (S  m- 1 ) 

r  tortuosity 

v\  volume  fraction  of  the  ionic  conductor 
(electrolyte  phase)  in  the  cathode 
vs  volume  fraction  of  the  electronic  conductor 
(catalyst  phase)  in  the  cathode 
01  potential  of  the  ionic  conductor 

(electrolyte  phase)  (V) 

0s  potential  of  the  electric  conductor  (catalyst  phase) 
(V) 

co  mass  fraction 


2.  The  model 

A  typical  free-breathing  cathode  together  with  a  layer  of 
electrolyte  is  illustrated  in  Fig.  1 .  It  is  common  to  DMFCs  or 
microPEM  fuel  cells  [14,15].  A  staggered  array  of  circular  holes 
is  distributed  on  the  top  surfaces  of  the  module.  They  are  open 
to  the  ambient.  A  current  collector  sits  against  the  remainder  of 
the  top  surfaces.  A  unit  cell  of  1 .5  mm  by  1 .5  mm  in  surface  area 
is  considered  (Fig.  2).  The  thickness  of  the  porous  cathode  and 
the  electrolyte  layer  is  the  same,  i.e.,  75  jam.  The  radius,  r,  of  the 
breathing  holes  varies  from  r  =  0.4  to  0.9  mm.  The  dimensions 


Fig.  1 .  Schematic  drawing  of  passive  cathode  with  a  perforated  current  collector. 

and  coordinate  system  of  the  computational  module  are  given  in 
Fig.  2. 

The  porous  cathode  contains  electrical  conductors  (catalyst), 
ionic  conductors  (the  electrolyte)  and  the  gas  mixtures.  The  oxy¬ 
gen  reduction  reaction  that  takes  place  in  the  porous  cathode  can 
be  expressed  by  the  following  equation: 

02  +  4H+  +  4e“  —*■  2H20  (1) 

The  breathing  holes  are  not  only  the  entrance  of  fresh  air  but 
also  the  exit  of  the  water  vapor. 

2.7.  Governing  equations 

No  reaction  takes  place  in  the  electrolyte  layer  that  connects 
the  two  electrodes  of  the  fuel  cell.  The  charge  conservation  in 
the  electrolyte  layer  is  given  by: 

V(-oiV0O  =  0  (2) 

where  o\  is  the  ionic  conductivity  of  the  electrolyte  and  <p\  is  the 
electrolyte  potential. 

In  the  porous  cathode,  the  current  that  passes  through  the 
cathode  is  composed  of  two  parts,  i.e., 

i  =  is  +  if  (3) 

is  and  ii  are  the  currents  flowing  through  the  catalyst  and  the 
electrolyte,  respectively.  Because  the  porous  cathode  is  elec¬ 
troneutral  everywhere,  no  charge  is  built  up.  Thus,  the  charge 
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conservation  becomes 

V  •  i  =  0  (4) 

That  is 

V  •  is  =  -  V  •  i  (5) 

These  two  current  components  interact  through  the  electrochem¬ 
ical  reactions.  What  is  transferred  from  the  electrolyte  to  the 
catalyst  is  denoted  by  the  (volumetric)  transfer  current  density, 
i.e., 

V  •  is  =  -  V  •  ij  =  -Sa  •  let  (6) 

where  Sa  is  the  surface  area-to-volume  ratio  of  the  porous  cath¬ 
ode  and  ict  is  the  local  transfer  current  density.  Thus,  by  using 
Ohm’s  law,  the  charge  conservations  in  the  catalyst  and  elec¬ 
trolyte  can  be,  respectively,  represented  by  the  following  equa¬ 
tions: 


V  •  (  ^s,effV</>s)  —  Sa  '  i'ct 

(7) 

V  •  (-<Tl,effV</>])  =  Sa  •  kt 

(8) 

where  0S  and  0 1  are  the  catalyst  potential  and  the  electrolyte 
potentials,  respectively.  crs?eff  and  a\  ^  are  the  effective  elec¬ 
tronic  and  ionic  conductivities  of  the  catalyst  and  electrolyte, 
respectively.  They  are  modeled  as: 


*ct*Sa^H2o 

IF 


(15) 


where  F  is  the  Faraday’s  constant.  The  effective  diffusivities  of 
the  oxygen  (Z)o2,eff)  and  water  vapor  (Dn2o,eff)  in  the  porous 
cathode  follow  the  Bruggemann  model  [18,19],  i.e., 


£>o2,eff  =  £zDo2  (16) 

Z^H20,eff  =  STDu2o  (17) 


The  velocity  vector  in  Eqs.  (11)  and  (12)  is  coupled  to  Darcy’s 
law  accordingly: 

u  =  —  —  S/p  (18) 

/x 

where  k  is  the  permeability,  /z  the  viscosity,  and p  is  the  pressure. 
The  continuity  equation  for  the  gas  flow  mixture  is  the  sum  of 
the  continuity  for  the  three  involved  species,  which  yields: 


V  (py^p)  =  ^(mo2  -  2Mh2o)  (19) 

The  relationship  between  the  electrochemistry  and  the  species 
transport  kinetics  on  the  cathode  surfaces  is  depicted  by  the 
Butler- Volmer  equation: 


<A> 


c02,ref 


exp 


4a  F 
~ RT 


°s,eff  —  CTs(l  Sc)  VS  p\,s 

(9) 

Tl.eff  =  Ol(l  -  ec)v\Pi,\ 

(10) 

where  vs  and  v\  are  the  volume  fraction  of  the  catalyst  and 
electrolyte  in  the  catalyst  layer,  respectively.  p^s  and  pn  are  the 
refer  to  the  catalyst  and  electrolyte  connection  of  the  catalyst 
layer,  respectively  [16,17].  It  is  noted  that  only  a  long-range 
connection  of  the  same  particles  stretch  from  the  catalyst  layer 
to  the  electrolyte  which  ensures  good  conductivity. 

The  species  conservations  for  oxygen,  water  vapor,  and  nitro¬ 
gen  in  the  gas  mixture  are,  respectively,  expressed  as 


f  ch2o 
VcH20,ref 


2 

exp 


4(1  —  a)F 
RT 


(20) 


where  a  is  the  symmetric  factor,  R  the  universal  gas  constant, 
and  T  is  the  temperature.  co2  and  ch2o  are  the  concentrations  of 
oxygen  and  water  vapor,  respectively. 


2.2.  Boundary  conditions 

The  boundary  conditions  for  the  ionic  current  balances  are 
insulating  at  all  boundaries,  except  for  the  bottom  of  the  elec- 


V  •  (pu&>02)  =  pZ)o2,e ffV&>o2  H - Z)° 2VM  +  Sq2  (11) 

poo  o2 

V  •  (pUOJu2o)  =  pF> H2o,effVu;o2  H - — Ah2oVM  +  Sh2o 

(12) 


&>N2  =  1  —  <J>02  —  ^h2o 


(13) 


where  p  is  the  density  of  the  gas  mixture,  u  the  velocity  vec¬ 
tor,  and  M  is  the  mole  masses  of  the  mixture.  In  Eqs.  (11)  and 
(12),  the  first  term  on  the  right-hand  side  represents  the  Fick 
diffusion,  while  the  second  term  is  the  correction  to  enforce 
the  Stefan-Maxwell  equations  on  the  multi-component  diffu¬ 
sive  system.  The  source  terms  of  the  mass  balances  for  oxygen 
and  water  vapor  are,  respectively,  given  by: 


*'ctSaMo2 
4 F 


(14) 


Table  1 

Electrochemical  and  flow  characteristic  data  used  in  the  present  model 


Expression  Symbols  Data 


i0  1.0  x  10-3  Am-2 

Sa  1.0  x  109  m-1 


Exchange  current  density 
Surface  to  volume  ratio 
Ambient  temperature 
Electric  conductivity  of  the  catalyst 
Ionic  conductivity  of  the  electrolyte 
Dynamic  viscosity  of  gas  mixture 
Porosity  of  the  cathode 
Permeability  of  the  cathode 
Reference  mole  fraction  of  oxygen 
Reference  mole  fraction  of  water 
Ambient  pressure 
N2  molecular  weight 
02  molecular  weight 
H20  molecular  weight 
Mass  fraction  of  02  at  inlet 
Mass  fraction  of  N2  at  inlet 
Mass  fraction  of  H20  at  inlet 


T 

353  K 

&  s 

1000  Sm"1 

(71 

5.0  Sm"1 

ll 

3.0  x  10-5  ms" 

£ 

0.4 

K 

1.0  x  10-9  m2 

c02,ref 

3.6641 

c02,ref 

3.6641 

Pin 

1.013  x  105  Pa 

Mn2 

0.028  kg  mol-1 

M0l 

0.032  kg  mol-1 

4^h2o 

0.018  kg  mol-1 

(7)02,in 

0.233 

(7)N2,in 

0.766 

fc»H20,in 

0.001 
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trolyte  layer.  The  condition  at  this  boundary  is  set  to  a  constant 
ionic  potential,  <p\  =  0.5  V.  The  boundaries  for  the  electrical  cur¬ 
rent  balance  are  those  that  encompass  the  cathode  sub-domain. 
The  boundary  conditions  are  all  insulating  except  for  the  cur¬ 
rent  collector  boundary  where  the  potential  is  set  to  zero, 
(ps=0. 

For  the  species  balances,  the  boundary  conditions  are  also 
those  that  encompass  the  cathode  and  are  all  insulating  except 
for  the  ambient-contact  boundary.  The  conditions  at  the  module 
inlets/outlets  (perforations)  are  given  by  the  composition  of  the 


fresh  air,  oo\  =  ca^o.  For  the  momentum  balance  through  Darcy’s 
law,  all  boundaries  are  no  slip,  except  for  the  module  inlet/outlet, 
where  the  pressure  is  set  top  =  /?atm.  The  input  data  to  the  model 
is  given  in  Table  1 . 

2.3.  Numerical  schemes 

The  governing  equations  are  numerically  solved  by  the  finite 
element  method  [19,20].  Computations  are  performed  on  8789 
unstructured  meshes.  A  grid-independence  test  is  carried  out 


(b)  r=0.75mm 


5  m&miui :  mi  si  $  c 
•  PPPfl fffj  ffm  ft  1 5 


(c)  r=0.9mm 


2*  ^  *  -t  .  ^  ^  %  '  v*  ‘  .  '  c - H  t: 

<  v  .A  ,  >  *  >  '  A  ~ 


Fig.  3.  Velocity  distributions  in  porous  cathode  for  three  different  perforations,  the  symbols  on  each  plot  represent  the  locations  of  the  maximum  velocity,  unit:  m  s 
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by  running  three  additional  coarser  and  finer  meshes.  A  typical 
simulation  requires  about  300  min  of  central  processing  unit  time 
on  a  Pentium  IV  2.8  GHz  PC. 

3.  Results  and  discussion 

Fig.  3  compares  the  gas-mixture  velocity  distributions  in 
the  porous  cathode  for  three  different  perforation  radiuses,  i.e., 
r  =  4.5mm,  7.5  mm  and  9.0  mm,  respectively.  All  plots  show 
outward  flows  of  the  gas  mixture  from  the  porous  cathode  to 
the  orifices.  It  seems  to  be  contrary  to  the  paradoxical  sense  that 
the  gas  flow  should  run  into  the  porous  cathode  from  the  ambi¬ 
ent.  Such  transport  phenomena  are  explained  as  follows.  In  the 
present  free-breathing  geometric  design,  the  orifices  act  as  not 
only  the  inlets  for  the  reactants  (O2)  but  also  the  outlets  for  the 
products  (H2O).  As  shown  in  Fig.  4(a),  the  oxygen  is  delivered 


diffusively  into  the  porous  cathode  through  the  orifices  and  then 
is  consumed  continuously  by  the  oxygen  reduction  reaction  (Eq. 
(1)).  Simultaneously,  water  vapor  is  produced  during  the  reaction 
and  should  be  removed  from  the  porous  cathode.  A  balance  of 
the  production  and  consumption  reveals  an  increase  in  the  bulk 
mass  of  the  gas  mixture.  Thus,  the  gas  mixture  in  the  porous  cath¬ 
ode  should  flow  outward  to  meet  the  continuity.  It  is  further  seen 
from  Fig.  3  that  significant  velocity  occurs  around  the  periphery 
of  the  orifice.  This  is  because  the  products  (water  vapor)  in  the 
reactive  layer  beneath  the  current  collector  accumulate  and  exit 
directly  from  the  periphery  of  the  orifice.  In  contrast,  the  gas- 
flow  velocity  near  the  orifice  center  is  rather  small  because  only 
the  products  beneath  the  orifice  are  removed  from  this  region. 
In  addition,  the  outlet  velocity  increases  with  decreasing  orifice 
radius.  A  smaller  orifice  is,  a  higher  is  the  outlet-flow  veloc¬ 
ity.  The  maximum  velocities  are  wmax  =  3.895  x  10-4ms-1, 


(b)  Mass  flux  of  oxygen 


(b)  Mass  flux  of  water  vapor 


Fig.  4.  Mass-flux  vectors  for  oxygen  and  water  vapor  in  the  porous  cathode  for  r  =  0.9  mm  the  symbols  on  each  plot  represent  the  locations  of  the  maximum  flux, 
unit:  Pasm-1. 
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2.187  x  10-4 ms-1, and  1.482  x  10~4  ms-1  for  the  perforation 
radius  of  r  =  0.45  mm,  0.75  mm,  and  0.9  mm,  respectively. 

Fig.  5  shows  a  comparison  of  the  iso-surfaces  for  the  oxygen 
mass  fraction  (coo2)  in  the  porous  cathode  for  three  different 
orifice  sizes.  Fig.  6  shows  the  corresponding  coo2  distributions 
on  the  bottom  surface  of  the  porous  cathode  (z  =  0.075  mm),  i.e., 
at  the  interface  of  the  electrolyte  layer  and  the  porous  cathode.  In 
these  two  figures,  the  total  overpotential  is  fixed  at  rjt  -  0.5  V.  It  is 
seen  from  these  plots  that  the  projection  region  from  the  orifice 
has  a  high  oxygen  mass  fraction.  It  decreases  as  the  oxygen 


gradually  diffuses  into  the  cathode  beneath  the  current  collector. 
Two  remote  corners  of  v  =  1 .5  mm,  y  =  0  and  v  =  0,  y  =  1 .5  mm 
have  a  local  minimum  value  of  coo2 ,  where  it  is  too  long  to 
breath  fresh  air.  As  for  the  effect  of  the  orifice  size,  the  oxygen 
mass  fraction  increases  with  increasing  the  orifice  radius  in  the 
entire  porous  cathode. 

Fig.  7  shows  the  effect  of  the  orifice  radius  on  the  overpoten¬ 
tial  distributions  on  the  surfaces  of  the  porous  cathode.  It  is  seen 
from  this  figure  that  the  surface  overpotential  decreases  with 
increasing  the  orifice  radius.  Near  the  orifices,  the  overpotential 


Fig.  5.  Distributions  of  iso-mass  fraction  surface  of  oxygen  for  r]t  =  0.5  V. 
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Fig.  6.  Distributions  of  oxygen  mass  fractions  on  the  plane  of  z  =  0.75  mm,  ijt  =  0.5  V. 
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Fig.  7.  Overpotential  distributions  on  the  module  surfaces  for  three  different  orifice  radiuses. 


is  small  due  to  the  lower  concentration  polarization.  In  contrast, 
the  remote  bottom  corners,  i.e.,  x=l.5  mm,  y  =  0,  z  =  0.75  mm; 
v  =  0,  pl.5mm,  z  —  0.75 mm,  have  significant  overpotential. 
Again,  this  is  attributed  to  the  lowest  oxygen  concentration  in 
these  regions.  Fig.  8  further  shows  the  overpotential  distribu¬ 
tions  along  the  distance  between  the  two  orifice  centers.  It  is 


seen  that  at  a  fixed  orifice  radius,  the  overpotential  has  the  low¬ 
est  value  at  the  center  of  the  orifice.  It  decreases  along  the  orifice 
radius  and  then  reaches  a  local  minimum  at  the  middle  location 
between  the  two  orifices  where  the  fresh  air  cannot  penetrate.  As 
for  the  effect  of  the  orifice  radius,  it  is  seen  from  this  figure  that 
the  overpotential  decreases  with  increasing  orifice  radius.  This 
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Fig.  8.  Overpotential  distribution  along  the  distance  between  two  centers  of  the 
orifices  for  z  =  0.075  mm. 

is  because  a  large  breathing  hole  can  reduce  the  concentration 
overpotential  in  the  porous  cathode. 

Fig.  9  shows  a  comparison  of  the  electrical  current  and  ionic 
current  distributions  across  the  vz  plane  at  y  =  0.  The  arrows 
indicate  the  direction  as  well  as  the  magnitude  of  the  currents. 
It  is  seen  that  in  the  porous  cathode  beneath  the  orifice  the 
electrical  current  is  directed  outward  from  the  orifice  center.  It 
increases  along  the  radial  direction  because  of  current  accumu¬ 
lation.  Near  the  radius  edge,  the  electrical  current  turns  upward 
to  the  region  that  subsequently  is  conducted  by  the  current  col¬ 
lector.  Beneath  the  current  collector,  the  electrical  current  goes 
upward  and  increases  gradually  from  the  bottom  of  the  porous 
cathode.  Again,  this  is  because  the  catalyst  receives  transfer 
current  from  the  electrolyte  during  the  electrochemical  reac- 
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Fig.  9.  Electrical  current  and  ionic  current  distributions  on  the  xz  plane  for  y  =  0. 


Fig.  10.  Effect  of  orifice  radius  on  the  electrochemical  performance. 


tion.  Thus,  the  electrical  current  accumulates  downstream.  As 
for  the  ionic  current,  it  is  quite  uniform  in  the  electrolyte  layer. 
In  the  porous  cathode,  by  contrast,  the  ionic  current  decreases 
upward  because  the  electrolyte  transfers  the  current  to  the  cat¬ 
alyst  by  the  electrochemical  reaction.  As  for  the  orifice  radius 
which  increases  from  r  —  0.45  mm  to  r  =  0.9  mm,  two  changes 
can  be  observed.  First,  a  large  breathing  hole  enhances  the  elec¬ 
trical  current  accumulation  beneath  the  orifice.  Secondly,  the 
ionic  current  is  enhanced  by  increasing  the  breathing  hole  from 
r  =  0.45  mm  to  r  =  0.9  mm. 

Under  a  fixed  total  overpotential,  the  electrochemical  per¬ 
formance  can  be  represented  by  the  current  density  that  passes 
through  the  fuel  cell.  Fig.  10  shows  the  effect  of  orifice  radius  on 
the  electrical  current  density  passing  through  the  computational 
module.  The  total  overpotential  is  fixed  at  0.5  V.  The  left  ordi¬ 
nate  is  the  total  current  drawn  from  the  cathode,  while  the  right 
ordinate  is  the  current  density.  For  a  common  basis,  the  current 
density  herein  is  averaged  over  the  entire  area  of  the  module 
(i.e.,  1.5  mm  by  1 .5  mm)  not  the  current-collector  area  only.  The 
curve  shown  in  Fig.  10  starts  with  an  increase  in  the  electrical  cur¬ 
rent  with  increasing  the  breathing-hole  radius  from  r  =  0.45  mm. 
It  reaches  a  local  maximum  at  about  r  —  0.65  mm,  and  then 
decreases  again  as  the  breathing-hole  radius  further  increases. 
As  the  radius  of  the  breathing  hole  increases,  two  effects  are 
present.  First,  a  large  opening  can  enhance  the  diffusion  in  the 
porous  cathode  that  reduces  the  concentration  polarization.  Sec¬ 
ond,  a  large  opening  reduces  the  area  of  the  current  collector  in 
contact  with  the  porous  cathode.  This  will  result  in  an  increase 
of  the  average  current  pathway  from  the  reacting  site  to  the  cur¬ 
rent  collector  that  increases  the  Ohmic  drop.  When  the  radius 
of  the  breathing  hole  is  around  r  =  0.65  mm,  these  two  compet¬ 
ing  effects  compromise  each  other  and  show  the  highest  current 
drawn  out  the  porous  cathode  under  a  fixed  total  overpotential. 

4.  Concluding  remarks 

The  effect  of  the  breathing-hole  size  on  the  electrochemical 
performance  of  a  free-breathing  cathode  of  a  DMFC  has  been 
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examined  in  detail  in  the  present  study.  A  three-dimensional 
model  that  couples  fluid-flow  fields  and  mass  transports  with 
the  electrochemistry  is  developed.  A  finite-element  based  CFD 
method  is  employed  to  solve  the  coupled  equations  with  a  proper 
account  of  the  electrochemical  kinetics.  Detailed  distributions  of 
electrochemical  species  characteristics  such  as  flow  velocities, 
species  mass  fraction,  species  flux  and  current  density  are  pre¬ 
sented.  A  larger  breathing  hole  increases  the  Ohmic  resistance 
penalty  because  the  electrical  current  has  a  longer  distance  to 
the  current  collector  from  the  reacting  site.  Simultaneously,  the 
area  in  contact  with  the  fresh  air  increases  with  increasing  r 
that  reduces  the  concentration  polarization.  The  two  competing 
effects  show  an  optimal  radius  of  about  r-  0.65  mm  that  pro¬ 
vides  a  highest  current  density  under  a  fixed  total  overpotantial 
across  the  computational  module. 
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